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ABSTRACT 



The "standard" expressions for total energy, linear momentum and also angular mo- 
mentum of asymptotically flat Bondi metrics at null infinity are also obtained from dif- 
ferential conservation laws on asymptotically flat backgrounds, derived from a quadratic 
Lagrangian density by methods currently used in classical field theory. It is thus a mat- 
ter of taste and commodity to use or not to use a reference spacetime in defining these 
globally conserved quantities. Backgrounds lead to Ncether conserved currents; the use of 
backgrounds is in line with classical views on conservation laws. Moreover, the conserved 
quantities are in principle explicitly related to the sources of gravity through Einstein's 
equations, while standard definitions are not. The relations depend, however, on a rule for 
mapping spacetimes on backgrounds. 
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1. Introduction 

In Science and Hypothesis, Poincare (1904) imagines astronomers "whose vision would be 
bounded by the solar system" because of thick clouds that hide the fixed stars. There is thus 
no fixed frame of coordinates, only relative distances and relative angles are measurable. 
In those circumstances, says Poincare, "we should be definitively led to conclude that the 
equations which define distances are of an order higher than the second. [...] The values 
of the distances at any given moment depend upon their initial values, on that of their 
first derivatives, and something else. What is that something else! If we do not want it to 
be merely one of the second derivatives, we have only the choice of hypotheses. Suppose, 
as is usually done that this something else is the absolute orientation of the universe in 
space, or the rapidity with which this orientation varies; this may be, it certainly is, the 
most convenient solution for the geometer. But it is not the most satisfactory for the 
philosopher, because this orientation does not exist." 

The conservation of energy, linear and angular-momentum, so useful in classical me- 
chanics and special relativity, are related to the homogeneity and isotropy properties at- 
tributed to an absolute "background" spacetime which does indeed not exist. In general 
relativity, the use of backgrounds is intrinsic to the definition of pseudo-tensor conserva- 
tion laws. Rosen (1958) and Cornish (1964) used backgrounds explicitly to calculate the 
energy-momentum in more appropriate coordinates than orthogonal ones. While back- 
grounds have thus proven to be useful, they nevertheless got little support and attention 
from the community of general relativists. On the contrary, great efforts have been spend 
to get rid of backgrounds, to avoid those "additional structure completely counter spirit 
of general relativity" (Wald 1984). True, efforts to avoid introducing background geome- 
tries generated interesting works notably by Penrose (1965), Geroch (1976) and other 
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mathematical physicists mentioned in Schmidt's (1993) review on Asymptotic flatness - 
a Critical Appraisal. As a consequence, we now have a rather well understood, coordi- 
nate independant, picture of asymptotic flatness, and we have also asymptotic coordinate 
independant expressions for globally conserved quantities in particular at null infinity. 

Conserved quantities at null infinity are given by integrals on spheres of infinite radius 
at fixed null time u = t — r = const. They represent total energy, linear momentum, 
angular-momentum and the initial position of the mass center (Synge 1964) of spacetimes 
with isolated sources of curvature at a moment of "time" u. The expressions have become 
"standard" according to Dray and Streubel (1984) [see also Dray (1985) and Shaw (1986)], 
and they are invariant under coordinate transformations of the Bondi-Metzner- Sachs, or 
BMS, group [Sachs (1962), Newman and Penrose (1966)]. Most standard expressions draw 
their strength - - beyond their esthetical appeal — from a few physically interpretable 
quantities, which have all been obtained before from the pseudo-tensors of Freud (1939) 
and of Landau and Lifshitz (1951). The quantities are 

(a) The Schwarzschild mass Ms c hw and more generally the dominant 1/r term in a multi- 
pole expansion of static solutions at spatial infinity (Geroch 1970) and at null infinity 
(Schmidt 1993), which was obtained from Einstein's pseudo-tensor [see Tolman 1934]. 

(b) The Bondi mass M(u) where u = t — r [see Bondi, Metzner and van der Burg (1962) 
see already Bondi (I960)] and Sachs's (1962a) linear momentum P(u) at null infinity* 
which have been calculated by M0ller (1972) using Freud's superpotential; 

(c) The Kerr angular- momentum L^err and more generally the dominant 1/r 2 "odd func- 
tion" factor [Misner Thorne and Wheeler (1973)] in a multipole expansion of stationary 
spacetimes at spatial or null infinity, which is obtainable from Landau and Lifshitz 

* What is actually defined in Sachs is dP/du. 
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superpotential for angular momentum [Papapetrou 1965]. 

The angular momentum at null infinity L(u), as given for instance by Dray and 
Streubel (1984), has not been derived from previously known soperpotential. L(u) is not 
a well defined physical quantity because of supertranslation freedom. However, the weak 
field approximation fits in well with results deduced from Landau-Lifshitz's superpotential 
(Creswell and Zimmerman 1986). The LL approximation is currently used in gravitational 
radiation calculations (Thorne 1980). 

Here we show that M(*u), P(u) as well as L(u) — and thus also Ms c hw and LKerr - 
can be calculated using a single superpotential derived from a quadratic Lagrangian with 
an asymptotically flat background at null infinity. The derivation insures automatically 
Poincare invariance, the absence of an anomalous factor of 2 in the M/L ratio and zero 
value for all conserved quantities when the spacetime identifies with the background it- 
self. As a result, dM/du and dP/du are the standard BMS invariant fluxes while dL/du, 
which is not BMS invariant, is nevertheless the same as the standard result. The position 
of the mass center, Z(u), associated with Lorentz rotations, and the corresponding flux 
dZ/du given here are not the same as the standard formula, at least in the non-linear 
approximation. 

Thus, the most important standard results, energy and angular momentum, can be 
deduced from Ncether conservation laws. One appealing aspect of Ncether conservation 
laws is that they are closer to classical intuition than the more abstract coordinate inde- 
pendant definitions at scri. Moreover, Ncether conserved quantities are directly related to 
the energy-momentum tensor of the matter through Einstein's equations by differential 
conservation laws. "Standard" definitions are not. 

Local differential conservation laws contain implicit definitions of local or quasi-local 
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conserved quantities which depend, on the local mapping of the spacetime on the back- 
ground. We give here no rule for local mappings. There are also several mapping- 
independent definitions of quasi-local energy in the literature, but they appear to be 
different from each other (Berqvist 1992). 

2. Superpotential and Noether Conservation Laws 

Let us briefly review the elements that lead to our superpotential (Katz 1996). Full details 
are given in a work by Katz, Bicak and Lynden-Bell (1996) — refered to below as KBL96 
- which has its roots in earlier work by Katz (1985) on flat backgrounds. KBL96 is on 
curved backgrounds and the superpotential obtained there is new even in the limit where 
the backgrounds become flat. 

(i) Lagrangian density for gravitational fields on a curved background 

Let g pu (x x ), A, fx, v, ... = 0, 1, 2, 3 be the metric of a spacetime M. with signature -2, and 
let g^ v {x x } be the metric of the background M. Both are tensors with respect to arbitrary 
coordinate transformations. Once we have chosen a mapping so that points P of M. map 
into points P of M., then we can use the convention that P and P shall always be given the 
same coordinates x x = x x . This convention implies that a coordinate transformation on 
M. inevitably induces a coordinate transformation with the same functions on M. . With 
this convention, such expressions as g^(x x ) —g^ v (x x ), which can be looked at as "pertur- 
bations" of the background, become true tensors. However, if the particular mapping has 
been left unspecified we are still free to change it. The form of the equations for "pertur- 
bations" of the background must inevitably contain a gauge invariance corresponding to 
this freedom. 

Let R x vpa and R vpa be the curvature tensors of M. and M.. These are related as 
follows [see Rosen (1940, 1963), see also Choquet-Bruhat (1984) for mathematical aspects 
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of background formalisms] : 

R\ pa = D P A}„ - D a A x up + A^AL - A^A^ + R\ pa . [2 - 1] 

Here D p are covariant derivatives with respect to g^, and A^ is the difference between 
Christoffel symbols in M. and M.: 

A x v = r*„ - f = \g Xp (D^g pu + D u g pp - D p g pv ). [2 - 2] 

Our quadratic Lagrangian density Lq for the gravitational field is then defined as 

C G = C-I , C = -^(R + d p k») , £=-^R , K = ^r- I 2 " 3 ] 

The mark " means multiplication by y/—g, never by y/—g, and a bar over a symbol signifies 
that g^v, Dp etc. are replaced by <7 M „, D p etc.. The vector density k p is given by 

& = ^=A,(-<hT) = g"A% - g p ° A%, [2 - 4] 

and its divergence cancels all second order derivatives of g^ v in R. t is the Lagrangian 
used by Rosen. C is C in which g pv has been replaced by g^. When g^ v = g pv , Lq is 
thus identically zero. The intention here is to obtain conservation laws in the background 
spacetime so that if g pv = g^, conserved vectors and superpotentials will be identically 
zero as in Minkowski space in special relativity. 

The following formula, deduced from [2-3] and [2-1], shows explicitly how Cq is 
quadratic in the first order derivatives of g^ or, equivalently, quadratic in A^ CT : 

C G = j-g^(A^A; a - A^Ay - - W)R^. [2 - 5] 

Notice that if R\ pa = and coordinates are such that f x u = 0, Cg is nothing else than 
the familiar "IT — TP' Lagrangian density [see for instance Landau and Lifshitz (1951)]. 
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(ii) Strong conservation laws and superpotential 
If 

Ax" = £ M AA [2 - 6] 

represents an infinitesimal one parameter displacement generated by £ M , the corresponding 
changes in tensors are given in terms of the Lie derivatives with respect to the vector field 
£ M , Ag^ v =£^g^ u AX, etc.. The Lie derivatives may be written in terms of ordinary partial 
derivatives <9 M , covariant derivatives with respect to g^ u , or covariant derivatives 
with respect to g^ v . Thus, 

= 9»\d u t x + g vX d^ x + i x d x g^ u [2 - 7a] 

= g^D„e + g vX D^ x + t, x D x g^ [2 - 76] 

= g»xD^ x + g„xD^ x . [2 - 7c] 

Consider now the Lie differential AC of C. With the variational principle in mind, we 
write AC =£^CAX in the form 

AC = ±G^Ag^ + «9 M >AA [2 - 8] 

where Einstein's tensor density, = — \g^ v R, is the variational derivative of C with 
respect to g Mi ,, is a vector density linear in £^ whose detailed form will not concern us 
here. The Lie derivative of a scalar density like C is just an ordinary divergence <9 M (££ M ), 
Thus 

6 = - d^C^) = 0. [2-9] 
Combining [2-9] with [2-8], we obtain 

6 = ±G^£zg^ + d»{A» - 1^). [2 - 10] 
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Bianchi identities imply D U G^ V = so that with [2- 7c], [2-10] can be written as the 
divergence of a vector density j^, say, 

6 = d^f = where f = \G^ V + B^. [2 - 11] 

Hence, £ "generates" a vector density j M that is identically conserved. It has been ob- 
tained without using Einstein's field equations; [2-11] is the kind of strong conservation 
law introduced by Bergmann (1949). We shall, of course, assume that Einstein's equations 
are satisfied, and replace -G£ by the energy-momentum of matter 

\G$ = T? [2 - 12] 

so that our strong conservation law [2-11] reads: 

d^ = d,(f^ u + B^)=0. [2-13] 

Equations [2-13] are, strictly speaking, not identities anymore. Given T^, [2-13] holds only 
for metrics that satisfy [2-12]. The vector density is linear in £ M and its derivatives up 
to order 2. 

Since j M as given by [2-11] is identically conserved whatever is g^ U: it must be the 
divergence of an antisymmetric tensor density that depends on the arbitrary g^s as well; 
thus 

f = d v f v , where f v = -f». [2-14] 

Indeed, is easy to find and is derived directly from C in Katz* (1985) [see also Chrusciel 
(1986), Sorkin (1988) and Katz and Ori (1990)]: 

= I £)[/*£!/] _|_ IflMfc"]. [2-15] 

* In the 1985 paper the background is assumed to be flat, but the derivation of does 
not depend on that assumption. 
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The terms ^_D^£^ will be recognised as | Komar's(1959) superpotential. In terms of D p 
derivatives, 

D p e = D p ^ + [2 - 16] 

and, using expression [2-4] for may be written in the form 

= 9 \mpd p ^ + g [ » P K[e + &g u]p K a - &A%g^, [2 - 17] 



Had we applied the identities [2-9] to C instead of £, we would have written everywhere 
~g pv instead of g^. We would have found strong, barred, conserved vector densities and 
barred superpotentials with the same £ M 's: 

> = f P C + B» = djji", [2 - 18] 

with 



= ku\M^] (k» = 0). [2 - 19] 

Strongly conserved vectors for Lq = C — C are obtained by subtracting barred vectors 
and superpotentials from unbarred ones; in this way we define relative vectors and in 
particular relative superpotentials J^ v — relative to the background space. Setting 

fn = jn _ JJ1 ; = _J^ = _ [ 2 _ 20] 

we have for the strongly conserved vector 1^ the following form: 

JM = JM + ( ?M = ^JM^ and dj" = 0, [2-21] 

which hold for any £^ and any mapping of Ai on A4; in [2-21] 

P = 6£?' + &»Md [p Z <T] . [2-22] 
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6%, a ppcr and are given explicitly in appendix for the interested reader. The relative 
superpotential density P v is now given by 

p v = ±(D^i^ - + | [ ^), [2 - 23] 

and can be also written in terms of g^, A^ a and £ M : 

K p» = i^D p ^ + g^ p A;[i x + iW"A% - e [M A^U PCT , [2 - n 

in which 

l* v = g pv - p 7 . [2 - 25] 

The tensors in [2-22] have a physical interpretation. On a flat background, in coordi- 
nates in which f * v = 0, [see appendix] 

#>=f^ + C [2-26] 

and reduces to Einstein's pseudo-tensor density. 9% appears therefore as the energy- 
momentum tensor of the gravitational field with respect to the background. The second 
tensor in [2-22], a^ pcr \ is quadratic in the metric perturbations just like t p . It is also 
bilinear in the perturbed metric components {g^ v — g^ v ) and their first order derivatives. 
fr»\p<r] resembles, in this respect, the helicity tensor density in electromagnetism. The 
factor of d[p£ a ] represents thus the helicity tensor density of the perturbations with respect 
to the background. 

It should be noted again that all the components of I p and of the superpotential 
J pu itself are identically zero if g^ v = g^; therefore strong conservation laws refer to 
"perturbations" only and not to the background. 
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(Hi) Ncether conservation laws 

We now consider what happens when arbitrary £ M 's are replaced by Killing vectors £^ of 
the background. defined in [2-22] , which contains the physics of the conservation laws 
[see KLB96], is not, in general, a conserved vector density since the identically conserved 
vector density is I M = + and thus 

d^P = -d^. [2-27] 

However, when £^ is a Killing vector of the background, then [see appendix] £ M = 
and J M (£) is conserved. 

Our has been derived in the same way as "Ncether's theorem" in classical field 
theory [see for instance Schweber, Bethe and Hoffmann (1956), or Bogoliubov and Shirkov 
(1959)]. Thus, by replacing £ M in strongly conserved currents by Killing vectors, £ M , of 
the background we obtain Ncether conserved vector densities in general relativity with 
mappings on curved backgrounds. 

>(£) = Kl v + ^ lpa] d [p L] = <V"(£) =0 [2 - 28] 

with J^ v (i) given by 

K j^(l) = l^D p e ] + 9 [ ^A;[e + £ [ V lp A^ - ^A^lgO*. [2 - 29] 

We can now integrate [2-28], on a part £ of a hypersurface S, which spans a two-surface 
<9£, and obtain integral conservation laws: 

c l V {t) = J m v + ^ W ^])^ = J j^(0d^„. [2 - 30] 

£ dT, 

V(£) depends only on the gravitational field and its first derivatives on dT, and on the 
mapping near the boundary. The relation with the matter tensor depends, however, on 
the mapping all the way down to the sources of gravity. 
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For weak fields on a flat background, the lowest order linear approximation of [2-30] 

is 

^V(0 = J Tf? = J J^{l)d^ v . [2 - 31] 

E <9E 

If the spacetime is asymptotically flat, intergrals over the whole hypersurface S ex- 
tending to infinity define globally conserved quantities. It is then appropriate to map the 
spacetime near infinity on a flat Minkowski space with its ten Killing vectors associated 
with spacetime translations, spatial rotations and Lorentz rotations. The ten Killing vec- 
tors give ten different expressions which can be interpreted respectively as the total 

— * — * 

energy E, the linear momentum vector P, the angular momentum vector L and the initial 
mass center position Z on S. 

We shall now calculate the conserved quantities for Bondi's asymptotic solution on a 
null hypersurface at infinity, using the right hand side of [2-30] . In what follows we intend 
to define all the quantities introduced. We shall not use cross-referencing for definitions 
which are often given with different symbols, factors of 2 or l/y/2 and other signs; most 
readers will appreciate this affort. 

3. Elements of the Asymptotic BMS Metric in Newman-Unti Coordinates 

(i) The Newman-Unti asymptotic solution 

In the coordinates used by Newman and Unti (1962) x x = (x° = it, x 1 = r, x 2 ,x 3 ), 
the metric of Bondi, Metzner and Sachs (1962) has the following form: 

T T 

ds = goodu + 2dudr + 2goLdudx + QKhdx dx , K, L = 2, 3. [3 — 1] 

The metric components are given by Newman and Unti in their formula (41). We shall 
make a few changes of notations because several indices which make sense in Newman and 
Unti are not useful here. Thus 
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* The indices of the metric where shifted from A = 1, 2, 3, 4 to A = 0, 1, 2, 3 (A or any 
other lower-case greek letter) so that x^u = x° = u, x 2 NU = x 1 = r, x% v = x 2 and 

4 3 

X NU ~ X ■ 

* We have denoted real parts with a prime like a', instead of Re(a) in NU, and imaginary 
parts with a "second" like a", instead of Im(a). Thus a = a' + ia" etc... 

* Complex conjugation is denoted here by an asterisk a*, because a bar over the symbol 
is reserved for the background. Thus U,nu = «*• 

* The ip° and a° of NU are here written ipj and a. Notice that ipj and a exist also in 
NU with a different meaning. 

* 0(r~ n ) of NU is here written O n . 

With these changes of notations, the metric components given by Newman and Unti 

become as follows: 

#oo = g 01 = 1 g 02 =g 03 = 

g- = -2P 2 VV* InP - ^ + iP'Mm-mV*") + 03 
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g +tg — j -\ 5 h C4 

# 23 = 3— + ^ + 6 

22 2P 2 4PV 6|a| 2 P 2 ^ 
2 = + ~ L 73— 

00 2P 2 4PV 6|a| 2 P 2 ^ 
2 =-72 ^ L 73— + 



[3-2] 



P is a function of x 2 , x 3 while i/>i, i/>2 and ex are complex scalar functions of it, x 2 , x 3 , 



defined in terms of the null tetrad components of the Weyl tensor. / is defined in terms 
of a and P 

f = 2P 4 V(a*/P 2 ), [3-3] 

and 

\/ = d 2 + td 3 . [3-4] 
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Note that we have expanded g 23 to the fifth order of 1/r. 
From [3-2] we have calculated the g^ v components: 

010 = 1 #11 = 912 = #13 = 

#00 = 2P 2 VV*lnP + 
1 



2^ 2 |P 2 V(^/P) + 1 IfP^ + aff 



902 + 1903 



IP 2 



923 = ~r 



a 



+ 



r 

1 <T (7 



r r J 

|Pa^ + 3|a| 2 r , " 
H 5 r C3 



p2 



+ 0, 



[3-5] 



p2 r p2 



+ o 2 



922 



933 = 



ra 



2P 2 P 2 2 P 2 



+ o 1 



r 2 ra' \a\ 2 ^ 
+ — - j^r + Oi. 



2 p2 p2 2P 2 

A useful quantity in our calculation is the density 



-9 = 



2P 2 



1 - 



[3-6] 



Further useful informations taken from Newman and Unti are 

(a) The w-derivative equations derived from the Bianchi identities, in the form given in 
their equations (40k, 1) or (42b, c). We shall most of the time use a dot on a symbol to 
denote a derivative of this function with respect to x° = u; thus ip = doifj. With this 
change of notation the formulas are 



V>i - PV^ 2 - 2ctV'3 = 



<A 2 + aa*-P 2 V(^1 =0 



[3-6] 
[3-7] 



where 



fa* 1 
\P 2 + P 



fa = -p3y J _ + __ V *y*p ) . 



[3-8] 

(b) The metric keeps the same form under coordinate transformations of the Bondi- 
Metzner-Sachs, or BMS, group (Sachs 1962b, Newman and Penrose 1966). The lead- 
ing terms in powers of 1/r of the transformation x x — > x x are given by NU in their 
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eq. (46): 



u 



r 



J(x L )u + K(x L ) +O l 
- T r + O 



[3 - 10] 



[3-9] 



x 



K 



Y K (x L ) + 1 . 



[3 - 11] 



J and K are arbitrary functions and the Y induce conformal transformations in 
(x K ) space, i.e. 



From this follows that P(x L ) can be fixed with an appropriate conformal transforma- 
tion [3-12], and r can be fixed by choosing a spherical boundary for the surface u = const. 
But u is only defined up to a supertranslation K{x L ). For a fixed P and a fixed r there 
are five independent scalar functions in the metric: and ipi", a' and a" and ip2 '• But 
they are defined up to a supertranslation K(x L ) and therefore there are actually four in- 
dependent initial quantities among the ten g^s on u = uq. The imaginary part of ip2, 
1P2" ', is not independent; It is defined in terms of P and a [NU (40g)]: 



The physical interpretation of these functions has been analyzed in details in a series of 
papers by Bondi, Metzner and Sachs [see especially Bondi van der Burg and Metzner (VII) 



d 2 Y 2 = ±d 3 Y s 



d 3 Y 2 = Td 2 Y 3 . 



[3 - 12] 




[3 - 13] 



1962] 
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4. Elements of the Asymptotic Background 

The asymptotic background is flat. In Minkowski coordinates X a = (X° = t, X k ) k, /, ... = 
1, 2, 3, the metric element 

ds 2 = riapdX^dX 13 = dt 2 - dX 2 . [4 - 1] 

An arrow designates spatial 3-vectors in Minkowski space. In X a coordinates, the Killing 
vector components of the ten translations are given by 

[4-2] 

and of rotations by 

efa/ 3] = (&7-Sw^ 7 - [4-3] 

In NU-coordinates x x = (u, r,x K ) 

ds 2 = g^dx a dx p = du 2 + 2dudr - — [(dx 2 ) 2 + (dx 3 ) 2 ] , [4-4] 

where 

r = a/a^ [4 - 5] 

and the metric of a sphere (u = r = r ) in conformal coordinates has the well known 
Riemann form for spaces of constant curvature [Eisenhart 1922]: 

P= I + I[( x 2 ) 2 + (x 3 ) 2 ] _ [4 _ 6] 

With P given by [4-6] we find that 

2P 2 VV*lnP = 1. [4-7] 

Eq. [4-7] somewhat simplifies g 11 given in [3-2]. Moreover, since P satisfies the following 
equation 

V*V*P = 0, [4 - 8] 
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ips defined in [3-8] becomes also simpler: 

/ has been defined in [3-3]. 

x 2 and x 3 are not uniquely defined nor are they uniquely related to the spherical 
coordinates {r,Q,<fi) m By definition X 1 = r sinO cos</>, X 2 = rsin6>sin</> and X s = 
rcos6>. We shall define x fl (X x ) as follows: 

x ° = u = t - r , x 1 = r [4-10] 
■P 



x- = V2—X 1 = v / 2cot(^) coscj) [4 - 11] 



„2 _ 

r 

x s = -V2^X 2 = -V2cot(^)sin0. [4-12] 



This choice makes the connection with the formalism of Newman and Penrose (1966) 
simple, as we shall see below. In terms of spherical coordinates that are sometimes useful 
in the calculations, P becomes 

p = Irr^- [4-13] 

2sin 2 (±6>) L J 

Let us introduce a unit vector in the radial direction in Minkowski space, denoted by 

— * 

e r = — = (sin6>cos<^>, sin6>sin</>, cos6>), [4 — 14] 

r 

and the complex 2-vector S on the sphere 

H = f + if. [4-15] 

In terms of e r and S, the ten Killing vector 4-components £ M or the 2 real + 1 complex 
components £ M = in {u,r,x K ) coordinates are as follows: 

(i) Time Translations 

£f 0) = ^ = (l, 0, 0), [4-16] 
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(ii) 3-Space Translations {tf } = A M : 

-e r , e r , — Ve r J , [4-17] 

(iii) 3-Space rotations {\e m kl ^ M] } = {rj^} = fp*: 

fp* = (0 , , 2iP 2 Ve r ) , [4 - 18] 

(iv) 3-Spacetime ("Lorentz") rotations = C M; 

C = (ue r , -(r + u)e r , - (l + ^) 2P 2 Ve r ) . [4 - 19] 

These Killing vectors satisfy the Killing equations 

D^ v + D v ^ = , ^=9^C [4-20] 

which are very useful in further calculations; using the background metric [4-4] , the Killing 
equations [4-20] can be written 



dl? = 


, do? = d 1 ? = -d ? 


[4- 


21] 


v(e°+e 1 ) = 




[4- 


22] 


ve° = 




[4- 


23] 


e _1 = 




[4- 


24] 


vs = 


0. 


[4- 


25] 



From equation [4-25] applied to fj^ defined in [4-18], and from equation [4-24] applied to 
C M defined in [4-19], we deduce, respectively, the following important identities 

V(P 2 Ve r ) =0 [4 - 26] 

P 2 V*Ve r = -e r . [4 - 27] 
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5. Globally Conserved Quantities for Bondi's Asymptotic Metric 

The hypersurface of integration is x° = const or u = uq, with boundary <9E at infinity: 
x 1 = r — > oo. On this boundary the conserved quantities V(£) defined in [2-30] may be 
written as surface integrals: 



no = ^ I J oi (0dx 2 dx 3 , [5-1] 



or in terms of the element of solid angle <ifi 



dx 2 dx i 



dtt = sin 0d6d(j) = ^ p2 [5 - 2] 



eq. [5-1] may be written as 



V(0= lim <f[P\j 01 m^. [5-3] 



kJ 01 {£) can be deduced from [2-29] or from [2-20], 



«j° 1 (e) = «5- 01 (e)-«5-° 1 (e)- [5-4] 

From [2-15], one obtains 

2«j 01 (0 = v rz 9(g la d l9aX e + d^ 1 - g la d a ?) + ^k 1 - f [5 - 5] 

which contains k° and A; 1 defined in [2-4]. The quantity j 01 is the same as [5-5] written in 
terms of g^ u rather than g^ u . Therefore, using g^ u defined in [4-4] and eq. [4-21], we have 



2«joi(e) = -2^=§d^ = -p2<9 £°; [5 - 6] 

remember that k° = k 1 = 0. 

We substitute the g^ and g^ u components according to [3-2] and [3-5] into [5-5], and 
calculate J 01 as given in [5-4]. We also use the Killing equations [4-21] and [4-24] to get 
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rid of We then place the resulting J 01 into [5-3] and take the limit r — > oo, to obtain 
the following expression for the integrant of V(^): 

f 



lim P z 2kJ ul (0 =<ra*doC - 

i — >oo 



2-02 + 2a&* + P Z V 



+ 



+ 



p2 



p 2 v 



i 



-o 1 



era - 



2P 2 



p2 



+ r 2pt [fE] '- 



[5 - 7a] 

[5 - 7b] 

[5 - 7c] 
[5 - 7d] 



The integrant of V(£) has been written as a sum of five terms. The term [5- 7c] gives 
no contribution to the integral [5-3] because — see [5-2] — P 2 disappears from [5- 7c] which 
becomes a pure divergence on a sphere, whose integral is zero. The term [5-7d] diverges 
as r — > oo. However, the factor of r integrated on the sphere is equal to the real part of 



1 ^dtt _ 1 
2~P2^"4^ ~ 27 



'E,dx 2 dx 3 , 



2P 4 



[5-8] 



which becomes, using [3-3] and [4-25] 



[5-9] 



If it is understood that one takes r — > oo after integration, the term [5-7d] does not 
contribute to [5-3]. Thus the integral of P 2 2kJ 01 contains only the factors of <9 £° and £° 
in [5- 7a], and the term with S in [5- 7b]. Further simplifications of [5-7] are still possible, 
but we want to keep it at present in this form. 

We shall now re- write the remaining terms of V(£) using the Newman and Penrose 
(1966) "edth" operators [see also Newman and Tod 1976], to enable comparisons with the 
formulas found in the literature of the 70's and 80's. More precisely, we shall denote by 
d the Newman and Penrose edth derivative times -^=. Thus with formula (3.9) of their 
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paper we define: 



drj = ^=d N p 66V = P^WiP'v) [5 - 10] 



where s is the spin weight of ?y, SW(?y) = s. The complex conjugate edth derivative 6* [see 
formula (3.17) in NP66] is defined by 

d*rj = P 1+S X7*(P- S 7]). [5-H] 

As far as our calculation goes, we must know the following spin weights: 
SW(ff) = 2 , SW(a*) = -2 , SW(dV) = 1 , SW(<9a*) = -1 

[5 - 12] 

SW(crff*) = , SW(e r ) = , SW(de r ) = 1 , SW(d*e r ) = -1. 
In terms of [5-10], [5-11] and [5-12], and with [3-3], we can re-write some of the quantities 
appearing in [5-7a,b] as follows: 

■i- = P 3 V(P _ V) = da* [5-13] 



P V 2^ = P V (, P 2P J = ^^P = 9 V [5 " M] 

= [5 - 15] 

PVe r = de r . [5 - 16] 

Eq. [4-26] and [4-27] can also be re-written as 

d 2 e r = [5 - 17] 

d*6e r = -e r [5 - 18] 

For integration by parts, the following property [equation (3.26) in NP66] will be particu- 
larly useful 

SW(A) + sw(S) = -1 j{6A)B d ^- = - j{dB)A^. [5-19] 
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With [5-13] to [5-16], the integral [5-3] of the integrant [5-7] may be written, taking 
account of [5-8] and [5-9] 



V(0 =\ i aa*d ^° 



4n 



4n 



[V>i + ada* - \d{aa*)\ 



2P 



4lY ' 



[5 - 20] 



Equation [5-20] for V(£) has a close ressemblance with the standard expressions build on 
a quite different basis. For comparison see for instance Winicour's (1980) formula given 
in [6-2] below. Equation [5-20] will now be simplified further. If we look at [4-17] and 
[4-19], we can see that the •u-components of the Killing vectors £° for space translations 
and Lorentz rotations are of the form £° = T{u)e r . For these £°'s, and with [5-17], the 
integral of the <9 2 cr*-term in [5-20] is zero; indeed 

-T{u) l[9Ve;ff = +F(u) Lda*de r ]'^ = -T{u) I \a* d 2 e r ]' ^ = 0. [5 - 21] 
J 4tv J 4tv J 4tv 

For the time translation the -u-component £^ = 1 (see [4-16]), and for this £° the integral 
of the <9 2 cr*-term in [5-20] is equally zero: 



1 

4n 



1 

8tt 



6 2 a*dQ = — & V J ^dx 2 dx 3 = 0. 



f_ 

2p 2 ' 



[5 - 22] 



The d a* -term does not contribute to space rotations for which ff is zero (see [4-18]) 



Therefore da* may be omitted from the integral [5-20], which reduces to 



V(0 =\ & aa*dot° — 



4tt 



11* + 



[il>! + ada* - hd(aa*)]=- 



2P 



dn 

4-k ' 



[5 - 23a] 
[5 - 236] 

[5 - 23c] 



Notice that S* = £ 2 — i£ 3 contributes only through their principal parts, with r — > oo. 
There is no in [5-23]. If we denote the principal parts of the 0,2,3 components of £ M 
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by £° = S( r=00 )) a = 0, 2, 3 we find, in the notations defined in [4-16] to [4-19] and in 
NU-coordinates that 



£( 0) = (i,o) 



[5 - 24a] 



A a = (-e r , 0) 



[5 - 246] 



rf = (0 , 2zPde r ) 



[5 - 24c] 




2Pde r ). 



[5 - 24d] 



The term d< 



>o£° = except for Lorentz rotations <f a , for which 



«9 C a = (e r ,0). 



[5 - 25] 



The £ a 's are the generators of the Poincare subgroup of the BMS transformations at null 
infinity (r = oo) [see Sachs 1962a, see also Newman and Unti 1966]. 

6. Detailed Comparison with Standard Results 

According to Dray and Streubel (1984), the expressions for angular momentum on the 
cross section S = {u = 0} given by Tamburino and Winicour (1966), Bramson (1975), 
Lind et al (1972), Prior (1977), Winicour (1968,1980) and Geroch and Winicour (1981) 
are all of the same form, though some of these authors differ by 'anomalous' factors of 
2 [see below]. For definiteness we shall compare [5-23] with the explicit expression for 
conservation laws L^(E + ) given in Winicour (1980) [his equation (2.16)]. The comparison 
is made easier with the following redefinitions of Winicour's notations, represented here 
with an index W: 



1968 paper , ip jw = 2V2ipj j = 1, 2 



d w = -V2d , a w = V2a , a w = 2d as u 



u 



[6-1] 



w 



V2~ 
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With [6-1], £ c (£+) can be writen 



in which 



and" 1 



L € (E+) = - fftih + ad* - 3*V)(£%)]'^ [6 - 2a] 

- 2 + ada* + ±6(vv*))(?ml)]'™ [6 - 2b] 

Z M = (1 , , , 0) [6 - 3a] 



"»m = ^(°. °. L -0; [ 6 - 3& ] 



£ M are the dominant parts of the asymptotic symmetry generators. Notice therefore that 
£ M m* in [6-2] is related to the conjugate of S defined in [4-15] as follows: 

= [6-4] 

(i) Energy and energy flux 

For time translations (see [5-24a]), we obtain from [5-23] an expression for the total energy 

E = V(Z m ) = -<f[ip 2 + aa*}'^-, [6-5] 

which is the expression given by Penrose (1964). Moreover, using [3-7] to eliminate ip2 we 
can calculate the flux dE/du of energy in terms of a: 

This is Bondi's (Bondi, van der Burg & Metzner 1962) mass loss formula. Both E and 
dE/du are also obtained from [6-2], using [5-22]. 



We are somewhat uncertain about the sign of m* , and were unable to decide. 
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(ii) Linear momentum and linear momentum flux 

For space translations (see [5-24b]), [5-23] provides an expression for the linear momentum 



P = T{\) = + j [V> 2 + a&*}' e 



r 47' 



[6-7] 



Which has also been proposed by Penrose (1964). Using again [3-7] and [4-7] the linear 
momentum flux is 

a a e r - — . [6 — 8] 



dP 

du 



4tt 



This dP/du is the BMS invariant flux first given by Sachs (1962b). Both P and dP/du 
follow also from [6-2], using [5-21]. 

(Hi) Angular momentum and angular momentum flux 

For spatial rotations in the background, defined by [5-24c], the corresponding conserved 
vector defined in the background is the angular momentum L 



L = V(ff) = 



(if)! + ada*)d*e r 



" dn 

4n ' 



[6-9] 



Notice that the term ^d(aa*)E* /2P, both in [5-22c] and in [6-2], drops out of the integral; 



indeed, using [5-19], we can write 



d{aa*) 



s* dn 

2P 4^ 



<9((7(j*)H<9*e r )^ 



i (f> (aa*)dd*e., 



dn 

' 4n 



[6 - 10] 



which is obviously zero because the integral is real. As a result, [6-2] gives twice (or perhaps 
minus twice) the value of L defined by [6-9] . The factor 2 is the well known anomaly pointed 
out by Penrose (1982) which makes that the Tamburino and Winicour (1966) as well as the 
Geroch and Winicour (1981) definitions do not agree with what one expect in the classical 
interpretation of the quantized weak gravitational field (Gupta 1952). The definitions of 
Streubel (1978), Dray and Streubel (1984), Dray (1985), Shaw (1986) and Penrose (1982) 
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do not have the anomalous factor 2. The angular momentum flux can be written in the 
following form using [3-6] and [4-9] 

^ = _ L(<j*d& - ada* + 2ada*)d*e r }"^. [6 - 11] 

du I 4tt 



(iv) The mass center initial position 

For Lorentz rotations in the background defined by [5-24d] and [5-25] , The conserved vector 



is 



Z = V(() = -uP+ <b {il) 1 + ad(j*)d*e 



4n : 



[6 - 12] 



because the term with <9o£ m [5-23a] cancels the term with d^aa*) in [5-23c]. The flux of 
Z is 



dz dp r „ 4 . A.^A^ v dQ 

= —u— <t> [(a da + ada )d e T 



du du 



An 



[6 - 13] 



Following [2-31], Z is the expression one expects in the weak field approximation. 
The conserved quantity Z w deduced from [6-2] is different: 



Z w = -uP + 2 



(V>i + ada* + \d\a\ 2 )d*e r 



' dtt 
4T' 



[6 - 14] 



remember the sign uncertainty of the integral. 
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7. Concluding Remarks 

V(^) defined in [2-30] is a coordinate independent expression. There are ten V{£), one for 
each of the ten Killing fields. For Poincare transformations in Minkowski space, the £'s 
and therefore also P(£)'s transform as vectors or tensors like in Special Relativity. Bondi's 
metric is, on the other hand, form invariant for supertranslations, [3-9] with J = 1 and 

~x K = x K 

u = u-K(x L ). [7-1] 

Such transformations induce changes in the values of the scalar functions a, ipi and ip2- 
In particular 

a(u,x K ) =a(u,x K ) -d 2 K. [7-2] 

The changes in dV/du are obtained with [7-2] from [6-6], [6-8], [6-11] and [6-13]. In 
particular, since 

we see that the fluxes of the energy and of the linear momuntum are unchanged by super- 
translations, except for a re-labelling of u into u; this is a well known result. 

We have thus found that the total energy, linear and also the total angular momentum 
at null infinity obtained from Ncether's theorem and a Lagrangian quadratic in first order 
derivatives, are the same as those considered as standard. The constructions of the stan- 
dard quantities are all, more or less, following the general principles outlined in Ashtekar 
and Winicour (1982). The use of Penrose's conformal spaces with the BMS symmetry 
group and the NP formalism is in some ways more appealing and more elegant than the 
background formalism used here. On the other hand, the mapping independant formula- 
tion has a serious drawback pointed out by Goldberg (1990): globally conserved quantities 
are unrelated to the source of gravity through Einstein's equations. While our background 
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metric formalism does not have this defect, it must be said that without a mapping rule the 
formalism remains incomplete, and we are stuck like everybody else with supertranslation 
ambiguities. Nevertheless, the fact that eq [2-30] relates V to the sources of gravitation is 
of potential value and fills, at least in principle, an important gap in the BMS invariant 
constructions on spheres at infinity. 
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Appendix 



See also KBL96 for detailed calculations. 



(i) 



[A -I] 



in which 



F v = g^ - g^, 



[A -2] 



and 



2< =g»° [(A x pX A^ + A» a Al - 2A" pX A* v ) - S^(A v pa A x x - Aj A Aj a ) 
+ ^ A (A^A^-A^A^). 



[A -3] 



(ii) 



is the antisymmetric part of a ppa and 



2k&^ = (g^g™ + 9^9 PU ~ <7 M V") A* A ~ {g vp g aX + 9 v °9 pX ~ 9 vX 9 pa )Kx- [A ~ 4] 



The two terms containing g pa do not contribute to a^ pcr \ 



(Hi) 



Ak^ ={Z^ p 9 pa + 9™Z° p - g^Z)A x x + (g^Z - 2g px Zl)A» pt 
+ F x d x Z + l pa (D»Z pa - 2D p Z t 



in which 



[AS] 



Z pa = £^g pa = D p ^ a + D^ p , Z = g pa Z pa and £ a = g a ^. [A -6] 
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